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ABSTRACT 



Context. A sample of 27 sources, catalogued as pre-main sequence stars by the Pico dos Dias Survey (PDS), is analyzed to investigate 
a possible contamination by post-AGB stars. The far-infrared excess, due to dust present in the circumstellar envelope, is typical for 
both categories: young stars and objects that have already left the main sequence and are suffering a severe mass-loss. 
Aims. The presence of two known post-AGB stars in our sample inspired us to seek for other very likely or possible post-AGB objects 
among PDS sources previously suggested to be Herbig Ae/Be stars, by revisiting the observational database of this sample. 
Methods. In a comparative study with well known post-AGBs, several characteristics were evaluated: (i) parameters related to the 
circumstellar emission; (ii) spatial distribution to verify the background contribution from dark clouds; (iii) spectral features, and (iv) 
optical and infrared colors. 

Results. These characteristics suggest that 7 objects of the studied sample are very likely post-AGBs, 5 are possible post-AGBs, 8 are 
unlikely post-AGBs, and the nature of 7 objects remains unclear. 

Key words, circumstellar matter — stars: post-AGB — stars: pre-main sequence 



1. Introduction 

In the study of large samples of stars having unconfirmed na- 
ture, the criteria used in selection of candidates often cause sam- 
ple contamination with objects of other nature than that of inter- 
est. This is the case of recurring confusion between pre-main se- 
quence (pre-MS) stars, and post-asymptotic giant branch (post- 
AGB) stars. In spite of their totally different evolution, both cat- 
egories of objects share common characteristics. The observed 
infrared (IR) excess, which originates from circumstellar dust, 
can be explained by re-emission of thermal radiation produced 
by the central source in both cases. 

Instigated by the presence of two confirmed post-AGBs in a 
sample of possible Herbig Ae/Be (pre-MS stars of intermediate 
mass), we decided to analyze in detail the objects having spectral 
energy distribution (SED) similar to those found in evolved stars. 
The sample was selected from the Pico dos Dias Survey (PDSfl 
by choosing only the PDS sources showing SED more luminous 
in the near-IR than in the optical, which is also a known charac- 
teristic of post-AGBs. 

Our goal is to distinguish, among the selected PDS sources, 
the very likely, the possible and the unlikely post-AGB objects, 
which could be included in "The Toruh catalogue of Galactic 
post-AGB and related objects'Q compiled by Szczerba et al. 
(2007). 

We first briefly describe the circumstellar characteristics of 
both, young stars and post-AGB objects in order to review simi- 
larities and differences that are relevant for the present work. 



1 PDS was a search for young stars based on infrared excess. 

2 http://www.ncac.torun.pl/postagb2 



1.1. Pre-main sequence stars 

Direct imaging is the most reliable way to study the geometry 
of the envelopes, and to establish input disk parameters for SED 
fitting models. However, sensitive imagery is constrained by sev- 
eral observational limits, being difficult to achieve for large sam- 
ples. 

In general, the circumstellar structure of pre-MS stars is 
traced through indirect means like spectroscopic data, for ex- 
ample, since the profile of spectral lines may be used to infer the 
physical conditions of line formation. 

Short-term spectral and polarimetric variability in Herbig 
Ae/Be (HAeBe) stars indicate, for instance, circumstellar non- 
homogeneities (Beskrovnaya et al. 1995), rotation, winds 
(Catala et al. 1999) or magnetic field (Alecian et al. 2008). The 
observed SEDs of HAeBes can also be used to explain their 
IR excesses, which are assumed to have originated from a disk 
and/or an envelope surrounding the central star. Different SED 
shapes have been used to classify the HAeBes according to the 
amount of IR-excess, which is one of the diagnostics of their 
evolutive phase in the pre-MS. 

The classification schemes of HAeBes (Hillenbrand et al. 
1992, Meeus et al. 2001) are based on the SED slope in the IR 
band, related to the amount of IR excess and the geometric dis- 
tribution of dust. However, these schemes do not consider more 
embedded objects, which correspond to the first phases of the 
young stellar objects evolution. A study of several other charac- 
teristics is required to check the pre-MS nature of the candidates, 
since these embedded objects have SED quite similar to those 
classified as post-AGB. 



1 



Vieira et al.: Post-AGB impostors in a sample of pre-main sequence stars 



1.2. Post-AGB objects 

Following Szczerba et al. (2007), in order to refer to the evolu- 
tionary stage after the AGB phase, we prefer to adopt the general 
term post- AGB instead of proto-planetary nebula, since a priori 
we have no information if our objects could become Planetary 
Nebulae. Post-AGB stars are luminous objects with initial mass 
between 0.8 and 8M Q . They completed their evolution in the 
asymptotic giant branch with a severe loss of mass (10 -7 to 10~ 4 
M /yr) (Winckel 2003). 

When the mass-loss process has removed all the envelope 
around the central core, the AGB star suffers a transition from 
highly embedded object to a new configuration of detached en- 
velope, being accompanied by changes in the SED shape, which 
acquires a double peak profile (Steffen et al. 1998). Van der Veen 
et al. (1989) classified a sample of post-AGBs, suggesting four 
classes according to the SED. Classes I to III have increasing 
SED slope from optical to far-IR wavelengths, while Class IV 
has double peak SED (maxima around near-IR and mid-IR). 

Ueta et al. (2000, 2007) studied the characteristics of de- 
tached dust shells and different SED of post-AGBs by means 
of J-K and K-[25] colors, which respectively describe the shape 
of the stellar spectrum in the near-IR, and the relation between 
the stellar (near-IR) and dust (mid-IR) peaks. According to 
these authors, post-AGBs can be separated in two kinds of mor- 
phologies: elliptical (SOLE) that seem to have an optically thin 
shell, with starlight passing through in all directions, and bipo- 
lar (DUPLEX) that have an optically thick circumstellar torus, 
where the starlight passes through only along the poles. Results 
from the HST survey of post-AGBs presented by Siodmiak et al. 
(2008) support this dichotomy in the morphology of the nebu- 
losities which was confirmed by Szczerba et al. (in preparation) 
for a large sample of post-AGBs. These authors found again 
clear differences in near-IR colors of SOLE and DUPLEX ob- 
jects, which are also correlated with the SED shape: post-AGB 
class IV objects (double peak) are SOLE type, while class II or 
III (single peak) are DUPLEX. 

Similar differences in SED shape are also found in pre-MS 
stars, indicating if the starlight is scattered in all directions (dou- 
ble peak SED) or not (single peak) (Hillenbrand et al. 1992, 
Malfait et al. 1998, Meeus et al. 2001). 

The main goal of the present work is to analyze a selected 
sample of possible HAeBes that is contaminated by the pres- 
ence of post-AGBs. As described in Sect. 2., the sample was 
extracted from the PDS catalogue by selecting the sources that 
show SED shape similar to post-AGBs. Seven characteristics, 
typical of evolved objects, have been used to distinguish them 
from the young stars in our sample. The analysis of these char- 
acteristics is presented as follows. In Sect. 3 we discuss the as- 
sociation with clouds in order to check (i) the effects of the inter- 
stellar medium in far-IR observed fluxes, and (ii) the occurrence 
of isolated objects. Spectral features and their relation with cir- 
cumstellar characteristics are discussed in Sect. 4. Optical and 
IR colors are evaluated and compared to known post-AGBs in 
Sect. 5. The main results are summarized in Sect. 6, according 
to the criteria used to identify the most probable evolved ob- 
jects, while concluding remarks and perspectives of future work 
are presented in Sect. 7. 



2. Sample selection based on IR-excess 

The Pico dos Dias Survey (PDS0) (Gregorio-Hetem et al. 1992, 
Torres et al. 1995, Torres 1999) was a search for young stars, 
based on the far-IR colors of T Tauri stars (TTs). Among the 
detection of more than 70 new TTs, PDS also provided a list of 
108 HAeBe candidates published by Vieira et al. (2003). 

Sartori et al. (2010) classified the HAeBe PDS candidates 
according to the SED slope (spectral index) and the contribu- 
tion of circumstellar emission to the total flux. They verified that 
84% of the studied PDS sources can be confidently considered 
as HAeBe, while the nature of several candidates remains uncer- 
tain. 

In this Section we describe the criteria adopted to select the 
objects studied in the present work, by using spectral index and 
the fraction of circumstellar luminosity as tracers of IR-excess. 

2. 1 . Spectral Index 

The PDS HAeBe candidates were separated by Sartori et al. 
(2010) into three groups, according to the SED slope between 
optical and mid-IR, measured by the spectral index /3\ = 
0.75 log(F l2 /F v ) - 1 (Torres 1999). We have selected 27 of the 
PDS sources of unclear nature, mainly those with high spectral 
index (J3\ > 0.7), which correspond to the most prominent IR ex- 
cesses. For objects showing single peak SED, the spectral index 
is related to the level of circumstellar extinction. Even if the SED 
is double-peaked, high values of the fi\ index also indicate high 
levels of circumstellar emission (at 12 fim), fairly larger than the 
optical flux (V band), independently of the stellar temperature. 

Among the selected sources there are two previously known 
post-AGBs: Hen 3-1475 (PDS465) and IRAS 19343+2926 
(PDS58 1), according to evidence reported in the literature (Riera 
et al. 1995, Rodrigues et al. 2003, Bowers & Knapp 1989), 
which motivated us to study their differences and similarities 
when compared to other objects of the sample. 

The list of the selected objects is presented in Table 1, giving 
their identification and the parameters used in the present work 
as spectral type (when available), B-V excess, IR colors (using 
2MASS, AKARI and IRAS data), and equivalent width of Ha line. 
Intrinsic polarization obtained by Rodrigues et al. (2009) is also 
given, when available. 

2.2. Fraction of Circumstellar Luminosity 

Different evolutive classes of pre-MS stars are defined according 
the observed SED, which is affected by several characteristics of 
circumstellar envelopes like: chemistry and size distribution of 
grains, and inclination of the system. 

The conspicuous grain features around 3.1 pm and 10 pm 
can be checked through ISO or Spitzer data, in order to better 
infer the nature of the circumstellar matter. However, due to the 
lack of near-IR spectral data for our whole samplfl disk mod- 
els cannot be constrained in the present work. For this reason, 
we decided to adopt a simple model (Gregorio-Hetem & Hetem 
2002) for the sole purpose of estimating the fraction of circum- 
stellar luminosity. Regardless of envelope or disk parameters, we 
are interested in determining the integrated observed flux, in or- 
der to estimate the contribution of circumstellar flux as a fraction 
of the total flux, defined by / Sc = (S, , a i-S s , ar )/S tota i. 

PDS was conducted at the Observatorio do Pico dos Dias, which is 
operated by the Laboratorio Nacional de Astrofisica/MCT, Brazil 
4 Only PDS141, 465, 518 and 581 were observed by ISO or Spitzer 



2 



Vieira et al.: Post-AGB impostors in a sample of pre-main sequence stars 

Table 1. List of studied stars 
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018 


05534254-1024006 


B7 


1.42 


2.17 


1.68 


7.41 


10.26 


2.04 


0.42 


0.10 


40 




027 


07193593-1739180 


B2 


1.41 


1.47 


1.32 


8.15 


9.78 


2.39 


0.84 


0.50 


88 




037 


10100032-5702073 


B2 


1.60 


1.86 


1.40 


8.28 


10.34 


2.55 


1.04 


1.00 


105 


2.76 


067 


13524285-6332492 


B 


1.56 


1.59 


1.49 


6.23 


8.56 


1.51 





-0.24 


90 


2.73 


141 


12531722-7707106 


? 




2.32 


1.82 


7.95 


11.06 


2.43 


0.84 


0.22 


40 




168 


04305028+2300088 


F0 


2.14 


1.74 


1.19 


6.51 


8.53 


1.87 


0.19 


-0.45 


8 


5.84 


174 


05065551-0321132 


B3 


1.01 


0.76 


0.65 


10.38 






1.17 


1.95 


65 




193 


05380931-0649166 


B9 


1.34 


0.97 


0.99 


7.20 


8.18 


1.86 


0.16 


0.52 


12 


2.25 


198 


05385862-0716457 


F0 


1.13 


1.03 


1.01 


6.95 


7.60 


2.79 


1.70 


1.07 


8.5 




204 


05501389+2352177 


Bl 


1.10 


1.24 


1.22 


9.68 


10.73 


3.31 


2.00 


0.96 


245 




207 


06071539+2957550 


B? 


1.20 


1.18 


1.21 


8.20 






0.62 


1.77 


5 




216 


06235631+1430280 


B2 


1.19 


1.32 


1.28 


8.14 


9.59 


1.65 


0.33 


2.03 


200 




257 


07414105-2000134 


A 


0.89 


1.20 


1.21 


7.81 






0.35 


1.77 


15 




290 


09261107-5242269 


A 


0.66 


0.35 


0.28 


10.9 






0.26 


2.33 


-7 


1.74 


353 


12222318-6317167 


B5 


0.99 


1.31 


1.17 


8.28 


9.68 


1.97 


0.45 


2.00 


200 


2.17 


371 


13473141-3639495 


09? 


1.40 


1.74 


1.31 


7.17 


9.18 


2.36 


0.58 


-0.1 


40 




394 


15351712-6159041 


F0 


0.44 


0.51 


0.52 


11.45 


11.17 


3.34 


1.51 


-0.7 


-2 


2.25 


406 


16050392-3945034 


A5 


0.40 


1.08 


1.05 


8.33 


9.19 


1.95 


0.79 


1.94 


13 


3.2 


431 


16545918-4321496 


AO 


0.52 


0.17 


0.16 


13.88 






2 


2.29 


8 


0.23 


465 


17451419-1756469 


B 


1.23 


1.36 


1.48 


8.37 


9.73 




1.51 


0.88 


110 


7.88 


477 


18003031-1647259 


Bl 


1.37 


1.43 


1.19 


7.42 


9.08 


2.22 


0.41 


-0.01 


120 


1.4 


518 


18273952-0349520 


OB 


2.30 


1.74 


1.35 


6.85 






0.37 


1.53 


600 


3.8 


520 


18300616+0042336 


F3 


1.26 


1.23 


1.09 


6.88 


8.13 


2.21 


0.70 


0.47 


33 


3.67 


530 


18413436+0808207 


A5 


0.47 


1.59 


1.66 


7.93 


9.92 


2.12 


0.52 


-0.07 


28 


11.06 


543 


18480066+0254170 


BO 


1.79 


0.48 


0.40 


8.91 


8.93 


3.82 


1.92 


0.41 


0.8 


1.14 


551 


18552297+0404353 


BO 


1.92 


1.53 


1.26 


8.16 


9.49 


2.15 


0.57 


0.57 


50 


11 


581 


19361890+2932500 


Bl 


0.89 


1.98 


1.72 


8.58 


10.82 




1.33 


0.74 


200 


12.22 



Columns description: (1) PDS name; (2) 2MASS identification; (3) spectral type; (4) B-V excess; (5, 6, 7, 8, 9,10, 11) infrared colors; (12) H„ 
equivalent width (negative values represent absorption lines); (13) intrinsic polarization. 

(Vieira et al. 2003), which makes uncertain the distinction be- 
tween a true association from a projection effect. Due to this 
restriction, we decided to evaluate the probability of association 
with clouds, where the objects apparently closer to dark clouds 
probably have young nature (age < 2 Myr), while stars located 
far from clouds are expected to be evolved objects. 

However, this cannot be used as a deterministic characteris- 
tic, since several examples of isolated pre-MS stars are known, 
as AB Aur, HD 163296, and HD100546, for example (e.g. van 
den Ancker 1999). Our main goals in investigating the spatial 
relation to dark clouds are twofold: (i) to verify the occurrence 
of isolated objects in our sample; and (ii) to estimate a possible 
far-IR background contamination. 

3.1. Distance to the edge of nearest clouds 

Aiming to infer the association with dark clouds, we have made 
use of the catalogue information compiled by Lynds (1962) and 
by Feitzinger & Stiiwe (1984). These complementary works al- 
lowed us to find the dark clouds closer to the objects of our sam- 
ple. For each dark cloud of these catalogues is assigned an opac- 
ity class number, which ranges from 1 to 6. We have adopted the 
criterion of choosing the closest dark cloud with opacity class 
3 at least. This choice restricts the selected clouds to those with 
relevant extinction level, although it must be kept in mind that 
these opacity classes were defined based on visual inspection of 
photographic plates. Given the selection of the closest clouds, 
we estimated their area and distance to each star of our sam- 
ple. These data have allowed us to roughly estimate the distance 
of our objects to the border of the selected clouds, defined by 
dedge = D - 0.5 x a 1/2 , where D stands for the projected dis- 
tance from the star to the cloud's center. By this way, a negative 
value for d e d ge means that a star is found "inside" the angular re- 
gion enclosed by the cloud of area a. A square area was adopted 



Figure 1 shows the synthetic reproduction of observed SEDs 
of our sample, estimated from the blackbody emission of three 
components: a central star (S star ), a flat passive disk (S^), sur- 
rounded by a spherical envelope (S e ). Different temperature laws 
are adopted: T^/T 

star ^ 

(r d )- 0J5 for the disk (Adams & Shu 

1986), and TJ 

T^ar (^e) for the envelope (Rowan-Robinson 
1986), as used by Epchtein, Le Bertre & Lepine (1990) to re- 
produce the IR data of carbon stars. The observational data are 
from PDS (optical photometry), 2MASS, MSX, IRAS, AKARI and 
ISO Catalogues. Three curves calculated by the adopted model, 
without fitting purpose, illustrate the expected variation on the 
fraction of circumstellar luminosity. In the case of PDS465, for 
example, f Sc varies from 0.90 to 0.94, which is the typical dis- 
persion found in our sample, leading us to adopt an error of 0.02 
in the estimation of /s c . 

All objects of our sample have [3\ > 0.7 and /s c > 0.7 that 
correspond to large amounts of circumstellar emission, suggest- 
ing two possible scenarios: they are embedded pre-MS stars or 
they possible are post-AGBs. 

We are aware that it is important to check if the data points 
around 100 pm are reflecting the cloud and not the immediate 
vicinity of the star, due to the large IRAS beam size in these 
bands, mainly for objects lacking of ISO spectra or AKARI data, 
for example. In the next section we analyze the background con- 
tribution from clouds, aiming to avoid a possible contamination 
of interstellar matter in the far-IR data. 



3. Association with dark clouds 

The star formation process is typically associated with clouds, 
where the gravitational collapse of pre-stellar cores gives rise 
to the embedded protostars. Thus, the association with these 
clouds may be an indication of young-star nature. Unfortunately, 
only rough distance measurements are available for our sample 
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0123 0123 0123 



log[A (/xm)] log[A (/im)] log[A (/mi)] 

Fig. 1. Observed SED of our sample showing log(/lF / | in [Watt trT 2 ] vs. log(/l) in [//m]. Filled squares represent optical photometry, 
2MASS, AKARI, MSX and IRAS data, while dots are used to plot ISO spectra. Solid curves indicate the variation of the calculated 
SED. 

in this first order calculation that does not take into account the 3.2. Reddening 
actual shape of the cloud, which usually presents a filamentary 

geometry. Nevertheless, d edge quantifies a possible association since the information obtained from catalogues of dark clouds is 
with the dark cloud. related to arbitrary opacities, we adopted two other methods to 

evaluate dark cloud effects on the objects of our sample. First we 
used an extinction map, which quantifies the cloud obscuration. 
In the second method we estimated the color excess in order to 
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Fig. 2. Position of some of the objects (displayed by the central circle) in the Galactic extinction map produced by Dobashi et al. 
(2005). The halftone grey color varies from to 3 mag, while the contours show Ay - 1 .0 mag level. 



obtain the total extinction, due to both interstellar (clouds) and 
circumstellar (disk/envelope) extinction. These two methods are 
described as follows. 



3.2.1 . Extinction Map 

We have adopted the Galactic map of visual extinction derived 
by Dobashi et al. (2005) using the classical star count method 
on the optical database of the Digitized Sky Survey I (DSS). 
The mapj provides the value of the extinction Ay for the region 
defined by \b\ < 40°. For each object, we have taken the Ay 
value and its mean dispersion in a neighboring region of 1 square 
degree, aiming to quantify the fluctuation of this extinction level 
in the interstellar medium. This quantity reveals the presence of 
a possible progenitor dark cloud behind the selected object. 

The obtained extinctions and respective dispersions are given 
in Table 2. For illustration, some specific regions of Dobashi et 
al.'s map are shown in Figure 2. It must be stressed that a pos- 
sible association of our objects with clouds may be due to pro- 
jection effects, where the dark cloud and the star are just on the 
same line of sight but not physically associated. However, since 
we do not know the distances of our objects, we can only speak 
in terms of probability of association. Objects having low values 
of Ay probably are isolated, but other methods are required to 
confirm this information. 



3.2.2. Color Excess 

The visual extinction may be also estimated by means of the 
color excess that is caused by the absorption and scattering of 
the emitted light, occurring not only in the ISM but also (and 
sometimes predominantly) in the circumstellar environment of 
our sample. Aiming to estimate this total extinction, we evalu- 
ated E(V— I) instead of E(B - V), since the color excess on V - I 
is less affected by the ultraviolet excess (Strom et al. 1975). The 
intrinsic colors were selected from Bessel et al. (1998) for dif- 
ferent sets of effective temperature and surface gravity. 

Based on the spectral type given in the PDS Catalogue, 
an estimative of effective temperatures was obtained by adopt- 
ing the empirical spectral calibration provided by de Jager 
& Nieuwennhuijzen (1987). Considering that the luminosity 



Available at http://darkclouds.u-gakugei.ac.jp 



classes are uncertain for most objects of our sample, the color 
excess was estimated by assuming two possible classes: main se- 
quence and supergiant that are adopted to represent HAeBes and 
post-AGBs respectively. Then, for each pair of luminosity class 
and T e ff we established a log g value extracted from Straizys & 
Kuriliene (1981), enabling the most convenient choice of intrin- 
sic colors. The E(V - I) values were converted to E(B - V) by 
the relation derived by Schultz & Wiemer (1975): E(B - V) = 
E(V - /)/(1.60 + 0.03). A normal reddening R v = 3.14 + 0.10, 
which is a mean value of the interstellar extinction for several di- 
rections of the Galaxy studied by Schultz & Wiemer (1975), was 
adopted to estimate the total extinction that includes both ISM 
and circumstellar effects, given by Ay to , - 3.14 x E(B - V). 
These are minimum values, since anomalous extinction (Ry >4, 
for example) may be found in dense interstellar clouds (e.g. 
Savage & Mathis 1979). We estimate that Ay tot can deviate by 
about 50%. 



3.3. Background contribution 

Far-infrared emission may also reveal the contamination by in- 
terstellar material on the observed SED. Aiming to evaluate the 
effects of background contamination, we have made use of IRAS 
images at 100 //m, which is dominated by the ISM cirrus emis- 
sion. For each IRAS image that contains one or more of the se- 
lected objects, we have chosen a region of 1 square degree in 
the neighborhood of the point sources, in order to estimate the 
background mean density flux and its standard deviation for each 
object. These background regions were selected avoiding the 
contamination of point sources in a radius of 1°. Table 2 gives 
the mean fluxes at 100/mi (<F>) obtained in the background re- 
gions. 

To quantify the excess of 100^m flux we compared the flux 
density of the point source (Fioo) with the mean background flux 
explained above, defining the fraction /fiqo = Fioo/(Fioo+<F>) 
as an indication of flux excess. If /fioo ~ the ISM dominates 
the local emission, while /fioo ~ 1 corresponds to an infrared 
excess mainly due to the circumstellar contribution. 

We also evaluated the fraction of "net" visual extinction, 
which roughly corresponds to the optical depth of the circum- 
stellar shell, by subtracting the interstellar obscuration (Ay) from 
the total extinction (Ay, ot ). In this case, we define /a v -{Ay tot - 
Ay)/Ay, ol , were the interstellar extinction is obtained from the 
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Fig. 3. Optical spectra of part of the sample, identified by the PDS number in the right side. The Ha feature is not shown in order to 
enhance details of other emission lines. Intensity scale is arbitrary. 



map of Dobashi et al. (2005) and the total extinction is calculated 
from the color excess estimated from magnitudes and spectral 
types given in the PDS catalogue (see Sect. 3.2.2). In this case, 
the effects of the adopted class of luminosity in the color excess 
estimation are negligible. 

Table 2 shows the estimated values of /a v and /fioo- Upper 
limits of /fioo were estimated for sources having bad quality of 
100 /mi data. Excepting PDS518, all the sources with well de- 
termined /fioo tend to show circumstellar extinction increasing 
with high levels of source emission. The two known post-AGB 
objects of our sample have /fioo > 0-5 an d /av > 0.5, meaning 
that circumstellar material has respectively prominent infrared 
flux and optical depth, significantly above the levels found in the 
corresponding background. We consider these high levels of /a v 
and /fioo characteristics of post-AGBs. 



4. Spectral Features 

Optical spectra, obtained by the PDS tearrfl are shown in Figure 
3 for the spectral regions containing [OI] (6302, 6365 A) and 



6 The original files were made available by S. Vieira, on the website 
http://www.fisica.ufmg.br/~svieira/TRANSF/ 



[SII] (6718, 6732 A) nebular lines, among others. The spectra 
of some sources are not displayed due to different reasons: they 
have no emission lines, excepting Ho- (PDS 174, 198, 290, 394, 
520, and 543); they are too noisy (PDS 168, 207, and 431); the 
6300-6500 A range was not covered (PDS 18, 27, 37, 67, and 
141), or the spectrum is not available (PDS257). 

Almost half of our sample shows the [OI] 6302 A emission 
line, but only PDS371 shows features similar to those found 
in evolved objects, like planetary nebulae (Pereira & Miranda 
2005). Considering the absence or unclear detection of these fea- 
tures for the whole sample, they were not used to diagnose the 
nature of our objects. 

The equivalent width of the Hq. line (Whq) of the objects 
in our sample was measured whenever possible and reported in 
Table 2. In the case of spectrum showing low signal-to-noise, we 
used the Wh„ published in the PDS catalogue. 

Most of the spectra show strong H„ emission line, but two of 
them appear in absorption (PDS290 and 394). Line variability is 
a possible explanation to the weak Hq,, usually seen in pre-MS 
stars. 

A comparison of the fraction of circumstellar luminosity 
(/sc) with the strength of H a line is shown in Figure 4. The dis- 
tribution of objects with Wh„ < 50 A spreads along the whole 
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Fig. 4. Equivalent width of H a line compared to /s c (top) 
and /av (bottom) for our stars, excepting PDS518, which has 
W(H Q ,)=600 A. Filled triangles are used to show the PDS sources 
with double peak in the SED, while open triangles represent sin- 
gle peak. 



range of /s c . On the other hand, objects showing strong emission 
line are mainly concentrated in the region of /s c > 0.87. As il- 
lustrated in Fig. 4, a similar distribution is found in the diagram 
of W Hq . in function of f Av , which is indicative of optical depth of 
the circumstellar shell, discussed in Sect. 3.3. 



5. Analysis of color-color diagrams 

To understand the nature of our objects it is necessary to carry 
out a multi-wavelength analysis. For this purpose, we have made 
use of two-color diagrams in several spectral ranges. The photo- 
metric optical data were obtained from the PDS catalogue. The 
measurements were taken in the Johnson-Cousins photometric 
system UBV(RI) C , with the high-speed multicolor photometer 
(FOTRAP) described by Jablonski et al. (1994). 

The catalogues 2MASS, IRAS, MSX, and AKARI (when avail- 
able) provided the data for the infrared colors analysis of our ob- 
jects. Optical and near-IR photometric data were corrected for 
extinction based on the A\r lot calculated in Sect. 3.2.2, adopting 
the A A /A V relations from Cardelli et al. (1989). 

5.1. Optical colors 

The spectral region represented by the B — V vs.U — B diagram 
is dominated by the photospheric radiation, which makes it very 
interesting to study the nature of the central source. The optical 
colors diagram is shown in Figure 5 displaying those of our ob- 
jects that have available UBV magnitudes. For each object, two 
different pairs of colors are plotted according to the extinction 
correction explained in Sect. 3.2.2 (by adopting two luminosity 
classes). Superimposed to this distribution are plotted the theo- 
retical curves for log g-2,3, and 4, calculated by Bessel et al. 
(1998). 

It is expected that HAeBes would have log g - 4, similar to 
main sequence stars, whereas the post-AGB would have log g < 
2, typical of supergiants. However, as shown in Figure 5 the the- 
oretical curves are clearly separated only for positive B - V val- 
ues, due to the presence of the Balmer discontinuity, which is 



B-\ 



Fig. 5. Optical colors of the studied sample, after deredden- 
ing using different luminosity classes: main sequence (filled 
squares) and supergiants (open squares). The intrinsic colors 
from Bessel et al. (1998) are displayed for different surface grav- 
ities. 



sensitive to surface gravity in this spectral range, being more 
significant for supergiants. In this case, only PDS406 seems to 
satisfy the gravity criterion. On the other hand, for B - V < 0, 
where part of our objects are located, the intrinsic colors from 
models for log g- 3 and 4 are similar. In this case, we consider 
U - B excess typical of evolved objects the values above the the- 
oretical lines, as indicated for PDS465 and 581, the post-AGBs 
of our sample, as well as PDS174, 216, 353, 477, and 518. 

5.2. Infrared colors 

The IR colors probe the thermal emission from grains that can 
be originated from dust at different temperatures: (i) cold grains 
found in circumstellar shells at temperatures lower than 300K, 
or (ii) warm dust (T ~ 1000K) produced by a recent episode of 
mass loss in the case of post-AGBs, or from a protostellar disk 
in pre-MS stars. 

5.2.1. Far-IR 

Diagrams of IR colors have been used to identify post-AGB can- 
didates according to their expected locus that is related to their 
evolution. Van der Veen & Habing (1988) suggested an evolutive 
sequence in a diagram of IRAS colors for post-AGB stars, de- 
fined by their mass-loss rate, in a scenario of short period Miras 
evolving into long period OH/IR stars. However, this hypothesis 
is based on observational results that are not firmly supported 
since the distances to OH/IR stars are uncertain (Whitelock et 
al. 1991). Other authors (e.g. Feast & Whitelock 1987) favor 
other scenario in which the Mira period is a function of its initial 
mass, where little changes in period or luminosity occur during 
the Mira evolution. This second scenario is supported by the re- 
sults, obtained by Whitelock et al. (1991), on the evolution of 
Miras and the origin of their Period-Luminosity relationship. 

Garcia-Lario et al. (1997), hereafter GL97, analyzed a sam- 
ple of post-AGBs comparing their position in the [12]-[25] 
vs. [25]-[60] diagram with several categories of objects, as for 
example: OH/IR variables (Sivagnaman 1989); pre-MS stars 
(Harris et al. 1988); and compact HII regions (Antonopoulos & 
Pottasch 1987). 
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[12]-[25] 

Fig. 6. Far-IR color-color diagram showing the regions for dif- 
ferent categories of objects (GL97). Filled triangles are used to 
show the PDS sources with double peak in the SED, while open 
triangles represent single peak. 



[9]-[18] 

Fig. 7. Mid-IR color-color diagram based on AKARI fluxes. 
Triangles indicate our sample and grey circles represent the very 
likely post-AGBs from the Toruh Catalog. Filled symbols are 
used for objects having double peak in the SED, while open 
symbols represent single peak. The distribution of young objects 
is shown by grey squares (YSOs detected by Spitzer) and grey 
crosses (HAeBe stars from HBC). 



Figure 6 shows the distribution of our sample in the IRAS 
color-color diagram, with respect to the regions occupied by dif- 
ferent classes of objects, as proposed by GL97. About 70% of 
our objects are located in the left part ([12]-[25] < 1) of this 
diagram, coinciding with the pre-MS box or above the HII re- 
gion box ([25]-[60] > 1.5). Our remaining objects, among them 
PDS465 and 581, are found approximately in the right side of 
the diagram ([12]-[25] > 1). The 12 pm flux density distribution 
for IRAS sources located in this part of the diagram was stud- 
ied by van der Veen et al. (1989, see region IV in their Figure 
1). Based on the spatial distribution compared to the number of 
sources as a function of the 12 pm flux density, they proposed a 
separation of objects in this region of the IRAS colours diagram. 
According to van der Veen et al. (1989) young stars, which are 
associated to clouds, are those sources with F12 < 2 Jy, while 
post-AGBs have F 12 > 2 Jy. 

5.2.2. Mid-IR 

The fluxes at 9, 18,65,90, 140, and 160//mfrorntheA£4fl/All 
Sky Survey (Ishihara et al. 2010) have been obtained for 78% of 
our sample. These data were used in the SED fitting discussed in 
Sect. 2.2. Following Ita et al. (2010) we verified the position of 
our sample in the J - [18] vs. [9] - [18] diagram, expressed in 
Vega magnitudes. In Figure 7, the AKARI colours of our sources 
are compared with the distribution of the very likely post-AGB 
stars from the Torun Catalogue (Szczerba et al. 2007). The posi- 
tion of young objects is also displayed in this diagram, by plot- 
ting the YSOs from the Spitzer survey of young stellar clusters 
(Gutermuth et al. 2009) and the HAeBe stars from the Herbig 
& Bell Catalog (HBC, 1988). In spite of the better quality of 
the AKARI data, when compared to IRAS colours, the overlap of 
different categories also occurs in Figure 7, not distinguishing 
young objects nor evolved stars. Most of our sample is located 
in this overlap region, but three sources (PDS 204, 394, and 543) 
coincide with the region of post-AGBs, appearing more clearly 
separated, outside the main locus of the HAeBe stars, for exam- 
ple. 



5.2.3. Near-IR 

Near-IR photometry can be used to check whether the main 
source of emission is photospheric, nebular or due to the pres- 
ence of a circumstellar envelope. Figure 8 shows the J—H vs. H— 
K diagram studied by GL97 to compare AGB and post-AGB 
stars with different categories of objects, like YSOs (region III); 
main sequence (region I/II); T Tauri and HAeBe stars (region 
III/IV), and planetary nebulae. GL97 verified that about 2/3 of 
the planetary nebulae are found in the "nebulae box" (region V), 
defined by Whitelock (1985). The very likely post-AGB stars 
classified by Szczerba et al. (2007) are also plotted to illustrate 
the distribution of evolved objects with single or double peak in 
their SED. 

Bessell & Brett (1988) studied the IR colours of Long- 
Period Variables (LPV), carbon stars and late-type supergiants. 
In Figure 8 we plot the schematic areas presented by Bessell & 
Brett (1988, Fig. A3), where carbon-rich stars are enclosed by 
dashed lines and oxygen-rich LPVs fall in the area defined by 
continuous line. The post-AGBs of our sample are located in re- 
gion III from GL97, beyond the right-hand-end schematic area 
from Bessell & Brett. The only object coinciding with the region 
of carbon-rich stars is PDS55 1 . 

It can be noted in Figure 8 that several of our stars are found 
in the same region as post-AGBs having single peak SED, while 
our objects having double peak SED appear in region I. As men- 
tioned in Sect. 1.2, Siodmiak et al. (2008) verified in near-IR 
diagrams a clear separation of nebulosity morphological classes 
SOLE and DUPLEX .The J - H vs. H - K diagram in Figure 8 
shows the regions occupied by DUPLEX (J - H >V) and SOLE 
(lower left side, H - K < 0.75). It is interesting to note that, in 
this diagram, the locus of the "stellar-like" objects dubbed by 
Siodmiak et al. (2008) (lower right side) coincides with the re- 
gion IV, expected for TTs and HAeBe stars (GL97). 

Comparisons of SED shapes are also interesting to be dis- 
cussed on basis of near- to mid-IR excess, which is related to the 
morphology of circumstellar structure (Ueta et al. 2000, 2007). 
Figure 9 presents the J - K vs. K - [25] diagram where stellar- 
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Fig. 8. Near-IR color-color diagram indicating the regions for 
different categories of objects defined by GL97 (separated by 
thin full lines), Siodmiak et al. (2008) (thin dashed lines), 
and Bessell & Brett (1988) (thick lines). Triangles are used to 
show our stars and grey circles represent very likely post-AGBs 
from the Toruh Catalog. Double peak SEDs are indicated by 
filled symbols (post-AGBs class IV). Grey crosses represent the 
HAeBe stars from HBC. 
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Fig. 9. Near- and mid-IR colors diagram indicating the regions 
defined by Siodmiak et al. (2008) (dashed lines). Filled triangles 
are used to show those of our stars with double peaked SEDs, 
while open triangles represent single peak SEDs. 

like post-AGBs are in the left part ( K - [25] < 8), while to 
the right are DUPLEXes (upper side) and SOLEs (lower part). 
Siodmiak et al. (2008) also correlate the SED shape with mor- 
phological groups: post-AGB class IV objects (double peak) are 
SOLE type, while classes I, II or III (single peak) are DUPLEX. 
As can be seen in Figure 9, the distribution of our double peak 
sources coincides with the SOLEs region. We consider these 
sources, as well as those located in the DUPLEX region, pos- 
sible being post-AGB stars. 

In spite of the good correlations between near-IR and SED 
shape, the discussion on elliptical or bipolar nebulosity types is 
only tentative since the morphology is unknown for most of our 
objects. One exception is PDS204, for which Perrin et al. (2009) 
used IR imagery and polarimetric data to confirm the presence 
of an edge-on disk surrounded by an extensive envelope pierced 



by bipolar outflow cavities. The results of the present work in- 
dicate that PDS204 has characteristics similar to post-AGBs, in 
particular its circumstellar morphology inferred from J-K X K- 
[25] colors, which is suggested to be DUPLEX (bipolar). 

6. Discussion and Summary of Results 

Aiming to analyse the similarities and differences between our 
sample of PDS sources and post-AGBs studied in other works, 
we have established several criteria to diagnose if the selected 
objects are pre-MS or not. The main result of this comparative 
study is the nomination of the likely post-AGB stars, indicated 
by notes in Table 2. For each given note, the adopted criteria 
and respective values that we consider typical characteristics of 
evolved objects are summarized as follows. 

(a) The possible association with clouds was investigated by es- 
timating the projected distance to the edge of the nearest cloud. 
According to the results presented in Sect. 3.1, note (a) indicates 
the objects having d e d ge > 0.5°. 

(b) The circumstellar extinction, defined by /a v , was estimated 
by subtracting interstellar extinction, obtained from Dobashi et 
al.'s Ay maps, of the total visual extinction, obtained from the 
E(B - V) calculation. The fraction of source emission at 100/im, 
defined by /fioo, was also evaluated in order to verify possible 
interstellar cirrus contamination. Note (b) indicates the objects 
showing / Av > 0.5 and / F ioo > 0.5, as described in Sect. 3.3. 

(c) Our sample was also analyzed according to the spectral fea- 
tures, in particular the equivalent width of H a line. In Figure 4 
the strength of H a is compared to the fractions of circumstellar 
emission /s c and /a,. Note (c) indicates all the objects showing 
Wn a > 50 A that have / Sc > 0.87 and / Av > 0.5. 

(d) In Figure 5 (optical color-color diagram), seven objects have 
U - B > -1 and B - V < -0.2 that we consider indicative of 
U - B excess typical of post-AGBs. As discussed in Sect. 5.1, 
no difference on gravities could be verified in the case of high 
temperatures, but one object with B - V > 0.1 seems to have low 
gravity, a possible indicator of evolved status. 

(e) As mentioned in Sect. 5.2.1, the IRAS color-color diagram 
(Figure 6) shows eight of our objects with [12] -[25] > 1, all of 
them having F12 > 2 Jy. This criterion has been suggested by 
van der Veen et al. (1989) to distinguish young stars, associated 
to clouds, from post-AGBs located in this part of the diagram 
that typically show F12 > 2 Jy. 

(f) and (g): The distribution of our sample in diagrams of near- 
IR colors was checked according the suggestion by GL97 that 
establishes different regions of these diagrams for each evolu- 
tive post-AGB phase. The objects located in regions I and III of 
Figure 8 are indicated by note (f). The diagram of near- and mid- 
IR colors shown in Figure 9 indicates that our objects having 
double peak SED are placed in the region of SOLE type post- 
AGBs and single peak sources coincide with DUPLEX objects, 
as suggested by Siodmiak et al. (2008). According the results 
discussed in Sect. 5.2.3, objects with K-[25]> 8 have note (g). 

7. Conclusions and Perspectives for future works 

The whole set of adopted criteria suggests that 26% (7/27 
sources) of the studied sample are very likely post-AGB stars, 
which achieved five or more notes. Four objects having three 
or four notes are considered possible post-AGBs. As discussed 
below, PDS371 is also included in this list indicating that 18% 
(5/27) of the sample possibly are post-AGBs. Among the other 
objects having less than three notes, 8/27 sources (30% of the 
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Columns description: (1) PDS name; (2) fraction of circumstellar luminosity; (3, 4) nearest cloud and respective distance to the edge; (5) visual 
extinction from Dobashi et al. (2005); (6) "net" visual extinction; (7,8) background mean flux and fraction of source emission at 100/um. 
Notes: (a) d edge >0.5°; (b) / Av > 0.5 and / F i o > 0.5; (c) high values of W Ua , /av and / F i o; (d) U-B excess; (e) [12]-[25] >1 and F l2 > 2 Jy; near-IR 
colors typical of post- AGB : (f) J-H and H-K, (g) J-K and K- [25] . 



sample) are unlikely post- AGB since they were confirmed in the 
literature as pre-MS stars. The nature of the remaining 7/27 ob- 
jects is unclear. Individual comments are given in Appendix A 
to present recent results from literature related to the nature of 
our stars. 

PDS465 and 581 are well known post-AGBs. The other very 
likely post-AGBs indicated by us are PDS 27, 174, 204, 216, 
and 353. Results from other works support these suggestions at 
least for two of them: PDS 27 (Suarez et al. 2006) and PDS 174 
(Sunada et al. 2007). 

The unlikely post-AGB objects: PDS 18, 37, 141, 168, 193, 
198, 406, and 518 are confirmed pre-MS stars (see Appendix A). 
Other sources with less than three notes are: PDS 67, 207, 257, 
477, 520, 530, and 55 1 , which were not previously studied in the 
literature. 

Other objects whose nature remains to be confirmed are 
PDS 290, 394, 43 1 , 543 that we suggest to be possible post-AGB 
objects. As discussed in Section 4, PDS371 is included in this 
group since it shows spectral features similar to proto-PN. The 
AKARI colours, discussed in Sect. 5.2.2, were not used to nom- 
inate the post-AGBs, due to the position of our sample, mainly 
concentrated in the region of overlap of young and evolved ob- 
jects in Figure 7. However, it is interesting to note that three 
sources appear out of this overlap region. Indeed, PDS 204 (very 
likely post-AGB), and PDS 394 and 543 (possible post-AGBs) 
have colours [9]-[18] > 3 mag, which coincide with the locus of 
most of the post-AGB class IV, confirming the classification that 
we suggest for these PDS sources. 

We are aware that additional observations are needed to 
definitively classify the candidates, as mentioned by Szczerba 
et al. (2007) in their evolutive catalogue of post-AGB objects. 
Several techniques and instruments operating in different spec- 



tral domains can be used to confirm the evolved nature of our 
candidates. Detailed spectroscopic analysis is required to obtain 
a better determination of spectral type and luminosity class that 
could be evaluated through the spectral lines sensitive to changes 
in gravity (Pereira & Miranda 2007). Spectral synthesis, for ex- 
ample, would provide abundances of C and O, whose ratio is 
used in the classification of post-AGB stars (van der Veen et al. 
1989). Abundances of Ba and Sr are also interesting to be deter- 
mined, since the enrichment of elements produced by s-process 
is expected to occur during the post-AGB phase (Parthasarathy 
2000). However, a detailed spectral analysis of absorption lines 
is constrained to objects with intermediary to low T e ff, which is 
not the case of our sources that do not show these spectral fea- 
tures. It is more indicated, in this case, to study the emission line 
profile of He I A 5876, or the strength of [N II] A5754 relative to 
[O I] A 6302, that were used by Pereira et al. (2008), for exam- 
ple, to indicate the nature of a proto-PN candidate. Evidence of a 
carbon-rich chemistry can be also obtained through mid-infrared 
spectroscopy from 10 to 36 fim, as the example of the study by 
Hrivnak et al. (2009). 

Direct imagery remains the most interesting tool to reveal 
the morphology of the nebulae. Our sources show several indi- 
cations of having a considerable amount of circumstellar ma- 
terial that could be detected in high resolution images. Ueta 
et al. (2007), for example, studied the post-AGB circumstellar 
structures via dust-scattered linearly polarized starlight, by us- 
ing imaging-polarimetric data obtained with NICMOS/HST. It 
would be very interesting to develop a SED fitting based on more 
reliable circumstellar parameters and a disk model more realistic 
than that adopted in the present work. We intend to apply a radia- 
tive transfer model specifically developed to study the physical 
conditions like mass loss and envelope density of the post-AGB. 
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Our results indicate new likely and possible post-AGBs that 
open promising and exciting perspectives to continue studying 
these targets. 
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Appendix A: Individual Comments 

The literature has been searched for results confirming the nature 
of the studied objects, which are found for 15/27 of the sample, 
all of them consistent with the results of the present work. In 
this section we separate these references in three parts: (i) previ- 
ously known young stars; (ii) very likely post-AGBs; (iii) possi- 
ble post-AGB or pre-MS. 

A. 1. Young stars 

PDS18 is included in the study of YSOs showing IR nebula, 
presented by Conneley et al. (2007). A reflection nebulae associ- 
ated to PDS18 is catalogued by Magakian (2003). Additionally, 
Conneley et al. (2008) also verified the presence of an arc-shaped 
nebula (~ 8 arcsec to the northeast) around this object. 

PDS141 (DK Cha), one of the most studied object of our 
sample, is a Class I pre-MS star. Results based on Spitzer data 
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have been presented by Porras et al. (2007), Alcala et al. (2008), 
Ybarra & Lada (2009), Evans et al. (2009), for example. The 
disk structure was studied by Liu et al. (2007). CO emission has 
been reported by Otrupcek et al. (2000). 

Spitzer data are also provided for PDS168 (Luhman et al. 
2006, Furlan et al. 2008), which is considered a Class I protostar 
(White & Hillenbrand 2004, Terebey et al. 2009). The protoplan- 
etary disk was studied from H2 emission (Bitner et al. 2008) and 
from millimeter data (Andrews & Williams 2005). 

PDS193 is in the 3pm spectroscopic survey of HAeBes, pre- 
sented by Acke & van den Ancker (2006). PDS406 was included 
in the study of Class protostars by Froebrich (2005) and ap- 
pears in the reflection nebulae catalogue from Magakian (2003). 

Sunada et al. (2007) included PDS198 in their study of H 2 
maser, but no maser emission was found for this object, proba- 
bly due to its intermediate-mass (F0 type star). Acke & van den 
Ancker (2006) consider PDS198 as HAeBe. 

PDS518 is also a well studied young B star (Acke et al. 
2008), for which 3pm spectroscopy of protoplanetary disks is 
presented by Goto et al. (2009); the circumstellar structure was 
also studied by Alonso-Albi et al. (2009). Conneley et al. (2008) 
included PDS518 in a study of multiplicity of embedded proto- 
stars. 



stars (Smolders et al. 2010), Acke et al. (2010) conclude that the 
chemistry of hydrocarbon molecules in the circumstellar envi- 
ronment is mainly affected by stellar radiation field regardless of 
the evolutionary status. 

A.3. Possible post-AGB or pre-MS 

PDS37 is in the list of post-AGB stars studied by Szczerba et 
al. (2007), but it is considered YSO in the work of Mottran et 
al. (2007) and in the Spitzer mid-IR survey by Carpenter et al. 
(2008). CO emission has been reported by Urquhart et al. (2007). 

PDS371 has been included in the photometry and spec- 
troscopy for luminous stars catalogue by Reed (2005). However, 
no additional information is available than those given by the 
PDS Catalogue (Vieira et al. 2003). 

CO emission has been reported by Brand et al. (1987) for 
PDS290 and PDS353. 

Pestalozzi et al. (2005) detected in PDS431 methanol maser 
emission, which usually is considered an indicator of massive 
star formation as well as observed in Class protostars. CO 
emission has been reported by Urquhart et al. (2007). 



A.2. Post-AGB 

PDS27 appears in the list of Suarez et al. (2006) who presented 
a Spectroscopic Atlas of post-AGB stars. CO emission has been 
reported by Urquhart et al. (2007). 

The emission of H2O maser (usually found in massive pro- 
tostars) was not detected in the survey by Sunada et al. (2007) in 
the case of PDS174 that is a result consistent with our suggestion 
of post-AGB candidate for this object. 

As stated before, PDS465 and PDS581 are confirmed post- 
AGB objects largely studied in the literature. The most recent 
results for PDS465 are presented by Lee et al. (2009), De Marco 
(2009), Luna et al. (2008), Raga et al. (2008), among others, and 
for PDS 581: Steffen et al. (2009), Montez et al. (2009), Dennis 
et al. (2008), Sanchez Contreras et al. (2008), for example. CO 
emission associated to PDS58 1 has been reported by Urquhart et 
al. (2008). Both objects are listed in the Toruh Catalogue, classi- 
fied by Szczerba et al (2007) as IRAS selected source (PDS465) 
or reflection nebulosity (PDS581). 

PDS 394 is studied as post-AGB candidate by Suarez et al. 
(2006) and Szczerba et al. (2007). 

PDS204 (MWC778) is considered a peculiar object in the 
H 2 maser study of post-AGB (Suarez et al. 2007, 2006). CO 
emission has been reported by Kawamura et al. (1998) and 
Urquhart et al. (2007). Herbig & Vacca (2008) suggested that 
MWC778 has a SED slope of pre-MS embedded object, which 
has spectral type F or G, due to the presence of metallic absorp- 
tion lines. From IR imagery and polarimetric data, Perrin et al. 
(2009) confirm the presence of an edge-on disk surrounded by an 
extensive envelope pierced by bipolar outflow cavities. They ar- 
gue that MWC778 has high bolometric luminosity inconsistent 
with F or G type star. The results of the present work indicate 
that PDS204 has characteristics similar to post-AGBs, in par- 
ticular its circumstellar morphology inferred from J-K x K-[25] 
colors, which are found in DUPLEX (bipolar) post-AGBs. 

PDS216 has been included in the study of hydrocarbon 
molecules in Herbig Ae stars (Acke et al. 2010) aiming to verify 
the correlation of effective temperature with the 7.8 pm feature, 
when comparing their sample with low- and high-mass stars. 
Considering that similar results are also found for S-type AGB 
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